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<n : 

After introducing the formalism of the general space and time fractional 
Schrodinger equation, we concentrate on the time fractional Schrodinger equation 

O ■ 

. and present new results via the elegant language of Fox's H-functions. We show that 

the general time dependent part of the wave function for the separable solutions 
of the time fractional Schrodinger equation is the Mittag-Leffler function with an 
imaginary argument by two different methods. After separating the Mittag-Leffler 

, function into its real and imaginary parts, in contrast to existing works, we show that 

the total probability is < 1 and decays with time. Introducing the effective potential 
approach, we also write the Mittag-Leffler function with an imaginary argument as 

: 

^ . the product of its purely decaying and purely oscillating parts. In the light of these, 

m ; 

0^ | we reconsider the simple box problem. 

cn '. 

Cn ' PACS numbers: 03.65. Ca, 02.50.Ey, 02.30. Gp, 03.65.Db 
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■ I. INTRODUCTION 

•l-H ' 

£ : 

The history of fractional calculus can be traced all the way back to Leibniz. In a letter 
to L'Hopital (1695), Leibniz mentions that he has an expression that looks like a fractional 
derivative [1, 2]. Later, Euler notices that using his gamma function fractional integrals 
may be possible. However, a systematic development of the subject does not come until 
nineteenth century, where Riemann, Liouville, Griinwald and Letnikov play key roles [1-6]. 
On the application side, after it became clear that anomalous diffusion can be investigated in 
terms of the fractional diffusion equation, the situation has been changing rapidly [1-6]. The 
first applications to quantum mechanics started with Laskin in 2000, where he considered 
the path integral formulation of quantum mechanics over Levy paths and showed that the 
corresponding equation of motion is the space fractional Schrodinger equation [7-17]. In 
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2004, Naber [18] discussed the time fractional Schrodinger equation [18-24] and obtained 
the time dependent part of the wave function for the separable solutions as the Mittag- 
Leffler function with an imaginary argument, which he wrote as a purely oscillatory term 
plus an integral that can not be evaluated analytically by standard techniques [18]. The 
main arguments and the conclusions of the paper were based on the critical assumption 
that this integral is uniformly decaying, hence could be neglected for sufficiently large times. 
This affected the main conclusions of the paper by giving probabilities greater than 1. Later, 
others have also followed this critical assumption [22-24] . 

We evaluate this integral exactly by using the elegant language of Fox's H— functions [25- 
27] and show that the decomposition of the Mittag-Leffler function as the sum of its purely 
oscillating and purely decaying parts is not possible. Since this point implies dramatic 
changes in the conclusions of the previous works [18, 22-24], we reconsider the problem of 
the time fractional Schrodinger equation. 

After introducing the time and space fractional Schrodinger equation in Section II, in 
Section III we concentrate on the time fractional Schrodinger equation and discuss its sep- 
arable solutions. We show that the time dependence is given as the Mittag-Leffler function 
with an imaginary argument by two different methods. After separating the Mittag-Leffler 
function into its real and imaginary parts, we show that the total probability is always < 1 
and decays with time. We also discuss the asymptotic forms of the wave function. In Section 
IV, we discuss the simple box problem, and in Section V, we introduce the effective poten- 
tial approach. Using the effective potential, in Section VI we show that the Mittag-Leffler 
function with an imaginary argument can be written, not as the sum, but as the product of 
its purely decaying and purely oscillating parts. In Section VII, we discuss our results. 



To write the fractional Schrodinger equation, we use the fact that the Schrodinger equa- 
tion is analytic in the lower half complex t— plane, and perform a Wick rotation, t — > —it, 
on the one dimensional Schrodinger equation, 



II. FRACTIONAL SCHRODINGER EQUATION 




h 2 d 2 ^ 

2m dx 2 



+ V(x)*(x,t), 



(1) 
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to write 



1 / d x - 



This is nothing but the Bloch equation [6] 

9* - <9 2 ^ 1 

where -D = - — is the quantum diffusion constant and V(x) is the potential. We can now 
write the time and space fractional version of the Bloch equation as 

$D?V(x,t) = \b a ,phPB?J>{x,t) - W(x)V(x,t), < a < 1, 1< (3 < 2, (4) 

ft ft 

where is the Caputo derivative and R% is the Riesz derivative (Appendix A). When 
there is no room for confusion, we also write fiD" = and i?f = V 13 = J^-. We have also 
introduced a new quantum diffusion constant, D a ,p-> with the appropriate units. Performing 
an inverse Wick rotation, t — > it, we obtain the most general fractional version of the 
Schrodinger equation as 

J^*0M) = jrDaj U^Yvfot) - l jV(xMx,t). (5) 



A. Space Fractional Schrodinger Equation 

When a — 1, Equation (5) becomes the space fractional Schrodinger equation investigated 
by Laskin [7]: 

J^OM) = {bi, P (n^y*(x,t) - ^v(xMx,t), (6) 



which he wrote as 

ih^*(x, t) = -Dp [hVf *(x, t) + V(x)*(x, t). (7) 

Laskin called [ftV]' 3 the quantum Riesz derivative and Dp in the above equation is the 
quantum diffusion constant for the space fractional Schrodinger equation, where as j3 — > 2, 
Dp ->■ l/2m. 

B. Time Fractional Schrodinger Equation 

For (5 = 2, we get 

— tf(x,f) = jD a , 2 [h—j *(x,t) - l -V(x)V(x,t), (8) 
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or 

Oqi 5^2 "ft 

OpVfrt) = i a D a , 2 h^(x,t) - l -V(x)*(x,t). (9) 

By defining a new quantum diffusion constant, D a = Da^h, we can write the time fractional 
Schrodinger equation as 

"ft 

— t) = t) - 0, < a < 1, (f 0) 

where -D a — > H/2m as a — > 1. Note that there is a certain amount of ambiguity in writing 
the time fractional version of the Schrodinger equation [18]. We will address this point again 
in Section VII. 



III. SEPARABLE SOLUTIONS OF THE TIME FRACTIONAL SCHRODINGER 

EQUATION 

We now reconsider the time fractional Schrodinger equation. Assuming a separable so- 
lution of the form fy(x,t) = X(x)T(t), we obtain the equations to be solved for X(x) and 
T(t), respectively, as 

dJ^- ~ ^X{x) = X n X(x), (11) 
^^=^ a XnT(t), 0<a<l, (12) 

where the spatial equation [Eq. (11)] has to be solved with the appropriate boundary 
conditions, and A n is the separation constant. At this point, the index n is redundant but 
we keep it for the cases where A is discrete. Since in the limit as a — y 1, D a — y h/2m and 
A n — y —E n /h, for physically interesting cases we will take A n < 0. 

A. Time Dependence — Method I 

Using Equation (A7) in Appendix A, we take the Laplace transform of Equation (12): 

= i a \ n £{T{t)} ) (13) 



d a T(t) 1 
dt a J 



to write 

s a f(s) - s Q - 1 T(0) = i a X n f (s), (14) 
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where T(s) = £ {T(t)} . This gives the Laplace transform of the time dependence of the 
wave function as 

c.a-1 



T(s) = T(0)- 



s a - i a \ n 

= T(0)- S -4 • 

v ' 1 - t a X n s- a 

Using geometric series, we can also write this as 

oo 

m=0 

oo 

= T(0) i am Ks~ 



-ma—l 



(15) 
(16) 

(17) 
(18) 



m=0 



which converges for \i a \ n s a \ < 1. The inverse Laplace transform of this can be found easily 
to yield the time dependent part of the wave function as 



°° ■am\mj.am 

T(t) = T(0) T , n s 
n r 1 + «™ 

= T(0) V ( " j 

m=0 v ' 

= T(0)E a (t a X n t a ), 



(19) 

(20) 
(21) 



where E a {i a \ n t a ) is the Mittag-Leffler function with an imaginary argument. 



B. Time Dependence — Method II 

To find the time dependence of the wave function, an alternate path can also be taken 
by directly evaluating the integral [18] 

1 P +ioc e^s^ds' 



T{t) = T(0) 



2m 



■y—ioo 



i a \, 



(22) 



The integrand has a pole at «i = due to the numerator, and a pole at s 2 = i\][ a due to 
the denominator, hence can be rewritten as (Bayin [6] and supplements of Ch. 14, [18]) 



T(t) = T(0) 



asm cot r 
a t io 



e xt x a 1 dx 



x 2a — 2a cos air x a + a 2 



(23) 



For \ n > 0, the above expression is valid for all < a < 1. For \ n < 0, the second 
pole, s 2 = |A n | 1//ct e *(*-A*+7r/2) ^ Yias both real and imaginary parts. Since the branch cut is 
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located along the negative real axis, we have to exclude the a = | , |, ^, • • • values, that 
is, a = 2/(5 + An), n = 0, 1, 2, . . ., so that the pole lies inside the contour. We now write 
Equation (23) as 



T(t) = T(0) 



-F a {a-t) 



where 



F a (a; t) 



a sin an 



^x^dx 



x 2a — 2a cos air x a + a 2 



(24) 



(25) 



In previous works, the above integral was assumed to be a decaying exponential, hence 
neglected in subsequent calculations [18, 22-24]. This resulted in probabilities greater than 
1. We now show that an exact treatment of this integral, which can not be evaluated an- 
alytically by standard techniques, proves otherwise. Using the elegant language of Fox's 
H- functions, details of which are presented in Appendices B and C, we evaluate the above 
integral exactly and show that it has an oscillatory part along with a Mittag-Leffler function 
as 



.,l/o. 



F Q (a;t) = ±- 



a 



-E a (\ n i a t a ). 



(26) 



Thus, the second method gives the time dependent part of the wave function as 



T(t) = T(0) 



iXn t 



a 



a 



+ E a {\ n et a ) 



(27) 



To be consistent with the robust result of Method I [Eq. (21)], we pick the minus sign, hence 
the time dependent part of the wave function is again obtained as 



T(t)=E a (\ n i a t a ), 



(28) 



where without loss of any generality we have set T(0) = 1. 



C. Probability and Asymptotic Behavior 

1. Total Probability: 

Since the general time dependence of the separable solutions of the time fractional 
Schrodinger equation is established as the Mittag-Leffler function with an imaginary ar- 
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gument [Eq. (28)], we can write the series 



00 (\ i a t a \ 

Til) = /•.',,! A,,/"/" ) = J2 r(l • - -v l A "'" 7 "l < ! 



i>=0 



X n i a t a 



+ 



r(l + a) r(l + 2a) r(l + 3a) 
which can be separated into its real and imaginary parts as 



Ea(\,n a ) = E«(t) + iE'Jt), < a < 1, 



= E 



^ [cos !ft] f< 



A^ [sin v -f-\ t av 



T(l + au) ' T(l + av) 

When ct = 1, naturally, E a (\ n i a t a ) becomes the Euler equation, that is, 

T a= i (t) = cos X n t + i sin X n t 



(29) 
(30) 

(31) 
(32) 



(33) 
(34) 



The eigenvalues, A n , come from the solution of the space part of the Schrodinger equation 
with the appropriate boundary conditions [Eq. (11)]. Note that we can also write Equation 
(32) as 

T{t) = E a (A„ cos 1 ^ (^) e) + iE a (A n sin 1 ^ (^) t«) , (35) 
which in terms of if— functions becomes (Appendix B, [25]) 



Tit) 



+i Hli (-a„W/» (^) r 



(0,1) 

(0,l),(0,a) 





(0,1) y 








(0,l),(0,a)/_ 



(36) 



An important consequence of this result is that for the time fractional Schrodinger equa- 
tion, the total probability, f*™ \^(x,t)\ dx, which is a function of time: 

»+oo 



+00 



\V(x,t)\dx = \T(t)f / \X(x)fdx 

J —CO 

Kit) I 2 + \ E ' tt (t) 1 2 
\E«(t)\ 2 + K(t)\ 2 , 



\X(x)\ dx 



(37) 



is not conserved. The normalization constant is fixed by the total probability at t — 0: 

"+00 



f 

J — c 



\^(x,0)\dx 



\X(x)\dx = 1. 



(38) 
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2. Small Time Behavior: 



For small times, we use the H— function representation of the Mittag-Lemer function: 

Ea{z) = H% (-4f U o, a) ) , (39) 
and write the time dependence as 

m = H% (-\^t%f ma) ) . (40) 
Using Equation (C3) we again obtain the series [Eq. (32)] 

cos(3cr7r/2) 



T(t) 



COS mi") 



cos(mr/2) 

1+ r(i + a) {Xnt ) + r(l + 2a) 



^ 2+ rW (A/)3 + 



+ i 



sin(mr/2) /x ^ i sin(mr) ^ +c ^ 2 : sin(3mr/2) 3 + 



-(AX) + 



-(AX) J + 



r(l + a) v """ y ' T(l + 2a) v " y r(l + 3a) 
which for small times gives the total probability as 



(41) 



/+oo 
-oo 



\^(x,t)\ 2 dx = \T(t)\ 2 



+oo 



\X(x)rdx 



x + 2cos(mr/2)(A ra r) + ^ 
aT(a) 



For A n < 0, this makes the total probability < 1 for small times. 



(42) 
(43) 



3. Large Time Behavior: 

For large times, we use the expansion in Equation (C7) to write 

-1 



v=0 



(1 - va)(\ n i a t a ) l+u 



\ n i a t a r(l- a)(\ n i a t a ) 2 r(l - 2a)(\ n i a t a ) 3 
We can also write this in terms of its real and imaginary parts as 

cos(a7r/2) cos(mr) cos(3a7r/2) 



T(t) 



+ i 



X n t a 
sin(o;7r/2) 

Xnt a 



r(l- a)(\ n t«y r(l-2a)(A n t«) 3 



+ 



+ 



sin(a!7r) 



sin(3a7r/2) 
r(l - a){X n t a ) 2 + T{l-2a){X n t a ) 



+ 



(44) 
(45) 



(46) 
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In agreement with Diethelm et. al. [28], for large times, the total probability decays as 

l T W| 2 «3^- (47) 

Since the wave function is zero on and outside the boundary, the energy does not leak 
out. Due to the presence of the fractional time derivative, the particle inside the infinite well 
is not free according to the ordinary Schrodinger equation. In Section V, we introduce the 
effective potential in terms of the ordinary Schrodinger equation, where the imaginary part 
of the complex effective potential describes the dissipation process. Hence, the energy is 
basically exchanged within the box, between the physical agent that generates the effective 
potential and the particle. 

IV. THE BOX PROBLEM 

We now consider a particle in an infinite potential well, where 

, < x < a 

V(x) = I . (48) 

I oo , elsewhere 

For the separable solutions of the time fractional schrodinger equation [Eq. (10)], we have 
already determined the time dependent part [Eq. (28)] of the wave function, hence we can 
write fy(x, t) = E a (\ n i a t a )X(x), where the spatial part comes from the solution of Equation 
(11) with the appropriate boundary conditions. For the particle in a box with impenetrable 
walls, we solve 

D a d ^^ = \ n X(x), tt(0,*) = tt(a,t) = 0, (49) 
which yields the eigenvalues and the eigenfunctions as 

A n = -D a (-) , X n (x) = C sin (— x) , n = 1, 2, . . . , D a > 0. (50) 

Since we are considering the time fractional Schrodinger equation, non locality is in terms of 
time, which means that the system has memory. For the separable solutions, the differential 
equation to be solved for the space part of the wave equation is an ordinary differential equa- 
tion, hence the boundary conditions used in Eq. (49) are the natural boundary conditions 
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for impenetrable walls. Using the normalization condition f \^ n (x, 0)| dx = 1, we write the 
complete wave function as 

* n (x,t) = J- [sin (—x)} E a {\ n i a t a ). (51) 

The total probability is time dependent, which for small times behaves as 

J \t, l (x,t)\ 2 dx=\E«(t)\ 2 + \Ei(t)\ 2 (52) 

^ i _2cos (a ,/2)(|A„|^) +0(t2tt)| (53) 

and for large times decays as 

lim / \ty n (x,t)\ 2 dx ~ -1 >0. (54) 

This is in contrast with the previous works [18, 22-24], where for large times the total 
probability is constant and greater than 1: 

lim [ \^ n (x,t)\ 2 dx = \ > 1. (55) 
t^oo J or 

A. The Energy Operator and the Box Problem 

For the time fractional Schrodinger equation, we take the new energy operator as 

ft f) a 

£ =_.?.£_ ( 56 ) 
i a dt a v 1 

where the Hamiltonian is 

H a = -D a H— + V(x). (57) 
With this operator, the energy values are real but time dependent: 

E = / V*(x,t)—V(x,t)dx 

i a J K ' dt a v ' 

= -h\ n \E a (X n t a t a )\ 2 J \X(x)\ 2 dx 

= -h\ n \E a (\ n t a t a )\ 2 (58) 
For the box problem, the energy eigenvalues are now obtained as 

E n = --J %(x,t)—V n (x,t)dx, (59) 

mt)\ 2 + K(t)\ 2 ] (60) 



hir 2 n 2 D a 



a 2 
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which for small times are given as 



E _ h<K 2 n 2 D a t 2cos(a7r/2)(|A w lt tt ; 



a 2 \ aT(a) 
For large times (a ^ 1), the system dissipates all of its energy as 

Jim En ~ h7[2n l D « o. (61) 

In the limit as a — > 1, D a — >■ H/2m, the energy eigenvalues approach to their usual values 
predicted by the Schrodinger equation: 

Note that the ordinary energy operator, yields complex eigenvalues, which only in the 

limit as t — > oo approach to real values that are larger than their usual values in Equation 
(62) by a factor of 1/V 2 . 



V. THE EFFECTIVE POTENTIAL 

To gain a better understanding of the effects of the fractional time derivative, we look 
for an effective potential in the (ordinary) Schrodinger equation that yields the same wave 
function. For simplicity we will set V(x) = in Equation (11). We now use the following 
useful formula [18]: 

^2 =? or [%Dtnt)\ + [ °™-"° ' (63) 

which is particularly helpful in isolating and exploring the effects of fractional time deriva- 
tives in terms of an effective potential. Using separable solutions of the time fractional 
Schrodinger equation [Eq. (12)], we can write 



dt 



^D]- a T{t) + 



(64) 



The quantity inside the square brackets is nothing but the Riemann-Liouville derivative [Eq. 
(A5)], thus 

^ = i«\ n [r L D}-*T{t)] . (65) 
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To define an effective potential for the separable solutions, we also write the (ordinary) 
Schrodinger equation as 

dT(t) . . h 2 d 2 X(x) rTi/ , /XTA/Xm/X , , 

X(;r) = "2^^^ T(t) + V eff-(t)X(x)T(t), (66) 

and substitute the wave function found from the solution of the time fractional Schrodinger 
equation [Eq. (65) and Eq. (11) with V(x) = 0]: 

i« 8 A» K^A 1 ^)] = -£-±X(x)T(t) + V e// .(t)X(s)T(t), (67) 

This yields the effective potential as 

^,(^.-«A„»!M + |1Al. ,68) 

Since the time dependence is given as T(t) = E a (\ n i a t a ) , we can also write 

V *"-®- 1 hXn E a (\ n i a t a ) + 2^D~ a (69) 

In the limit as a — > 1, L> Q — >■ , the Mittag-Lemer function becomes the exponen- 

2m 

tial function, T(t) = e lXnt , hence the effective potential vanishes as it should. With this 
effective potential, the Schrodinger equation yields the same wave function as the time frac- 
tional Schrodinger equation and satisfies the same boundary conditions. In other words, the 
wave function found from the time fractional Schrodinger equation, satisfies the Schrodinger 
equation with the above effective potential. 

Using the if— function [25] representation of the Mittag-Leffler function (Appendix B): 

E a (X T rt a ) = H$ (-Xni a tX'Mo,a)) ( 7 °) 

and Equation (B7) for the Riemann-Liouville derivative of the if— function, we can also 
write the following closed expression for the effective potential: 

,,,„ M , r/-l,2 (\ .■a / a|(0,a) ) (0,l) \ 
T/ (A {t 1+a h\ n ) H 2,3{ t |(o,l),(o,a),(l-a,a)J H 2 \ n 

£ffX) = ~^~ a ^L AjV |(o,D \ + 2^D a (71) 

11 1,2 y A n L t 1(0,1), (0,a) J 

Using the symmetries of the if— function [25], we can further simplify the above equation 
as 

T/ ^ (^ + ^AQ^ 1 :2 1 (((-^) i/ ^) a i;;;;; i(1 _ aia) ) ^ An 

ff wf_ AW .|(W) ^ + 2mfV ^ 

-"1,2 ^ ra 1 l(0,l),(0,a) ) 
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or as 

V ^- {t) ~ ti- E a (\ n i a t a ) + 2^' (73) 

where 

E a , p (z) = T ( a Z k + /3 y a ' G C ' Re (°) > °' Re ^) > 0, ^ G C, (74) 



fc=0 



is the generalized Mittag-Leffler function, and E^iCz) = -E'a(z) [25]. 



A. Hamiltonian 

Since the effective potential is time dependent and complex, the energies corresponding 
d 

to the operator ih— are also time dependent and complex, that is, the ordinary Hamiltonian 
at 

operator: 

H = -££ + v -"-® < 75 > 

is not hermitian. However, complex potentials in quantum mechanics have found interesting 
applications in the study of dissipative systems [29-32]. 

Let us now write the effective potential for the time fractional Schrodinger equation 
explicitly. Using Equation (73) and the expression 

°'° ( " ' ) " h r(i + + * £ r(i + to) ' (76) 

we can write the effective potential in terms of its real and imaginary parts as 

Ve f f.(t) = V e R ff ,(t)+tV e I ff Xt), (77) 



where 



v r , v = ^ 2 A n _ fiA w sin(a7r/2)t a 1 
e//A J 2m,D a r(a) 

-^(f^- ra+W ^ 1 ^- (78) 



and 



ft 2 A n /iA n cos(air/2)t 



a-l 



//A ; 2////A, ~ r(a) 



i ftA 2 sWrv^ ( 1 I cos(a7r) 2cos 2 (mr/2) . 2a-1 

+ ftA n 8in(Q7r) ^ (1 - Q)r(Q) + - r(Q)r(1 - q) j t + . (79) 

The imaginary part of the complex effective potential plays an important role in describing 
the decay of the energy of the stationary states of the time-fractional Schrodinger equation. 
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VI. THE MITTAG-LEFFLER FUNCTION, E a (\ n i a t a ), AND IT'S PURELY 
DECAYING AND OSCILLATING PARTS 

Solving the Schrodinger equation [Eq. (66)] with the above effective potential, we now 
obtain the following expression for the time dependent part of the separable solutions of the 
time fractional Schrodinger equation: 



T(t) = 



exp ( ~ h I Vl ff .(t')dt' 



exp 



h 



h 2 x„t 



2mD, 



-- ! V e R f ,{t>)dt> 

a J 



(80) 



where A n is a separation constant. This is naturally equivalent to the solution found before 
[Eq. (28)], that is, T{t) = E a (\ n i a t a ). Note that when a = 1, the effective potential 
vanishes, hence the time dependence reduces to 

T(t) = exp — ] , E n 



v h J'~" 2mDi (81) 
In other words, the time dependent part of the separable solution of the time fractional 
Schrodinger equation, E a (\ n i a t a ) , can be written, not as the sum, but as the product of its 
purely oscillating and purely exponentially decaying/growing parts. Despite the fact that 
the effective potential corresponding to the time fractional Schrodinger equation is complex, 
the advantage of the effective potential approach is that one can use the mathematical 
structure of the Schrodinger theory. 
Finally, using the expansion 

K//(t) ^^_5^« r _ 1 + 0(t ^ 

'h\ n cos(«7r/2) ! 



+ 1 



r(a) 



the small time behavior of T(t), to lowest order, becomes 

% ( h\ n sin(a7r/2) 



A n cos(a7r/2) 
exp | 1 



exp 



T(t) ~ . ., ,. , 

\ al (a) 

For the box problem, where A„ = — D a (^ E ) 2 , 
decays exponentially, which for small times becomes 



(82) 



(83) 



n 



h \ aT(a) 

1,2, ... , D a > 0, the total probability 



2 - 2c ° s(a7r/2) |A„|^ + 0(t-). 



\T{t)\'~l- 



aY{a) 



(84) 



This is in accordance with our previous result [Eq. (43)] obtained from T(t) = E a (\ n i a t a ). 



15 



VII. CONCLUSIONS 

The fact that anomalous diffusion can be studied by fractional calculus has attracted re- 
searchers from many different branches of science and engineering into this intriguing branch 
of mathematics [1-6]. Since 60's, successful examples are given in random walk, anomalous 
diffusion, economics and finance. Introduction of fractional calculus into a certain branch of 
science is usually initiated by replacing certain derivatives in the evolution or the transport 
equations with their fractional counterparts. A general feature of these applications is that 
replacing a time derivative with its fractional counterpart basically introduces memory ef- 
fects into the system and makes the process non Markovian, while a replacement of a space 
derivative introduces global or non local effects. However, it is important to note that not 
all such effects can be covered by fractional generalizations of the basic equations [20, 21]. 

In this paper, we first established the time dependent part of the separable solutions of 
the time fractional Schrodinger equation as the Mittag-Lefner function with an imaginary 
argument by two different methods. This showed that the total probability is < 1 and 
decays with time. We also introduced the effective potential approach, where with the 
effective potential, the (ordinary) Schrodinger equation yields the same wave function, with 
the same boundary conditions, as that the time fractional Schrodinger equation. 

d 

For the energy, the standard operator used in Schrodinger theory, E = ih—, yields 

at 

complex values since the corresponding Hamiltonian is complex and non hermitian. On 

h d a 

the other hand, the operator E = — — , which follows naturally from the time fractional 

i a dt a 

Schrodinger equation, yields real but time dependent energies [Eq. (58)]. In the light of 
these, we also discussed the particle in a box problem, which is essentially the prototype of 
a device or a detector with internal degrees of freedom. With the operator E = — — t^, 
the energy eigenvalues are real and time dependent, and decay with the imaginary part of 
the complex effective potential [Eq. (79)] as 

K = ^^exp^|K / // X^'j. (85) 

Since the wave function is zero on and outside the boundary, the energy is exchanged within 
the box between the source of the effective potential and the particle. The imaginary part 
of the complex potential describes the dissipative processes involved. 

Complex potentials have been very useful in the study of dissipative processes in both 
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classical and quantum physics. Sinha et. al. [29] have extended the factorization technique 
of Kuru and Negro [30] to study complex potentials in classical systems that are analogues 
of non hermitian quantum mechanical systems. In the context of quantum device modeling, 
a general non hermitian Hamiltonian operator has been investigated by Ferry et. al. [31]. 
They have also discussed the general behavior of complex potentials with applications to 
ballistic quantum dots and its implications for trajectories and histories in the dots. Barraf 
[32] have used uniform complex potentials to model particle capture from an incident beam 
by quantum wells. We have pointed out that there is a certain ambiguity in writing the time 
fractional Schrodinger equation. For separable solutions, the time fractional Schrodinger 
equation that Naber [18] and Dong and Xu [22] used gives the time dependence of the wave 
function as E a (X n (—i) a t a ). However, for the total probability [Eq. (37)] and the energy [Eq. 
(58)], both conventions yield the same time dependence since 

\E a (\ n (-i) a t a )\ = \E a (\ n i a t a )\ , < a < 1. (86) 



Appendix A: Basic Definitions 



Riemann-Liouville Definition of Differintegral: 

The basic definition of fractional derivative and integral, that is, differintegral, is the 
Riemann-Liouville (R-L) definition: 

For q < 0, the R-L fractional integral is evaluated by using the formula 

(Al) 



r d«f 




_[d(t-a)]i_ 


r(-q) l 



L- J [t- If]-*- 1 f (fidlf, q<0. 



For fractional derivatives, q > 0, the above integral is divergent, hence the R-L formula 
is modified as [6] 

't 

[ 

i 

where the integer n must be chosen as the smallest integer satisfying (q — n) < 0. 



d q f 


d n 




[[d(t-a)]*_ 


~ dt" 


J(n-q)J a 



[t-t'\- {q ~ n) ~ l f{t')dt' 



q > 0, n > q, 



(A2) 



For < q < 1 and a = 0, the Riemann-Liouville fractional derivative becomes 



dm) 

dti 



1 



R-L 



d r l f(t')dr n 



T(l - q) dx J (t - t')i 



(A3) 



Caputo Fractional Derivative: 
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In 60's Caputo introduced a new definition of fractional derivative [1-6]: 



<Pf(t) 
dti 



a 



* fdf(r) 



T(l-q)J V dr ) (t - r) 



dr 



-, < q < 1, 



(A4) 



which was used by him to model dissipation effects in linear viscosity. The two derivatives 
are related by 



W(ty 








c 


dti 



r«/(o) 



, < q < 1. 



(A5) 



Laplace transforms of the Riemann-Liouville and the Caputo derivative are given as 
(Supplements of Ch. 14 of [6]) 

n-l 

, n - 1 < q < n, (A6) 



k=0 

£{°mm=s q m-Il s9 - k - 1 dkf{t) 



t=0 



fe=0 



, n — 1 < q < n, 



(A7) 



t=o 



-, and wherever there is need for distinction we use the abbrevi- 



where a D?/(t) - ^ _ a)]q , 
ation 'R-L' or £ C. Since the Caputo derivative allows us to impose boundary conditions in 
terms of ordinary derivatives, it has found widespread use. 
Riesz Derivative: 

In writing the space fractional diffusion equation or the space fractional Schrodinger 
equation, we use the Riesz derivative, which is defined with respect to its Fourier transform 
(Supplements of Ch. 14 of [6]) 



F{B*f(t)} = -\u>\ q g(u>), 0<q<2, 



as 



»+oo 



where g(u) is the Fourier transform of f(t). Note that 



(A8) 



(A9) 



(A10) 



Appendix B: Fox's H- Function 



In 1961 Fox introduced the if— function, which gives a general way of expressing a wide 
range of functions encountered in applied mathematics. H— function provides an elegant 
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and an efficient formalism to handle problems in fractional calculus. Fox's if— function is 
a generalization of the Meijer's G— function and is defined with respect to a Mellin-Barnes 
type integral [25-27]: 



h(s)z- s ds, (B2) 



1 

2tH 



c 



where 

1 J UU + i T (a 3 + Ajs) Y\U + i T (l ~ b 3 - B j8 ) ' 1 6) 

m,n,p,q are positive integers satisfying < n < p, 1 < m < q, and empty products are 

taken as unity. Also, Aj, j = 1, . . . , p, and Bj, j = 1, . . . , q, are real positive numbers, and 

a,, j — 1, . . . ,p, and bj, j — 1, . . . , q, are in general complex numbers satisfying 

Aj(b h + v) ^ B h ( aj - X - 1) for v, X = 0,1,...; h=l,...,m, j = l,...,n. (B4) 

The contour C is such that the poles of T(bj + Bjs), j = 1, . . . , m, are separated from the 
poles of T(l — cij — Ajs), j = 1, ... such that the poles of Y{bj + Bjs) lie to the left 
of C, while the poles of T(l — aj — Ajs) are to the right of C. The poles of the integrand 
are assumed to be simple. The if— function is analytic for every \z\ ^ when /i > 0, and 
analytic for < \z\ < 1//3 when /i — 0, where /i and /3 are defined as 

i=i j=i 

' n >/ n^/ (B6) 

j"=i i=i 

Fox's if— function is very useful in the study of stochastic processes and in solving 
fractional diffusion equations [33, 34]. For example, the following useful formula for the 
Riemann-Liouville fractional derivative of the if— function [33]: 



o iJ 



„a,Trm,n ( ( \b I (°j > A j)\ 1 _ a-f3 rj-m.n+l / / n 6 1 ( "'^("i^) \ (Ti7\ 
Z n p,q \} CZ > \{b j ,B j ))\- Z n p+l,q+l \\ CZ ) \{b j ,B j ),{p-a,b)) ' ^ ' J 



where a, 6 > and a + bmm(bj/Bj) > —1, 1 < j < m, can be used to find solutions to 
the fractional diffusion equation by tuning the indices to appropriate values. Similarly, the 
Laplace transform of the if— function can be obtained by using the formula [25] 

£ {^Ktu } = s ' p Kf (-"iS::t)) > (B8) 
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where the inverse transform is given as 

n-l / „-pTrm,n ( erl {o, v ,A v )\ \ _ „p-l Tjm,n ( -er I («p A0,0v)\ 
X Y U P,q \ aS \(b q ,B q )J / ~ X n p+l,q\ aX \(b q ,B q ) )> 

where 

p,a,s E C, Re(s) > 0, a > 

and 



Re(p) + <7 max 

KKn 



1 Re(aj 



7Tt 



> 0, |arga| < — , 9 = a — a. 
2 



(B9) 

(BIO) 
(Bll) 



Appendix C: H-Function in Computable Form 

Given an if— function, we can compute, plot, and also study its asymptotic forms [25-27] 
by using the following series expressions: 
I) If the poles of 

m 

Ylntj+sBj) (ci) 

are simple, that is, if 



B h (bj + A) ^ Bj(b h + is), j ^ h, j,h = l,...,m, A, v = 0, 1, 2, . . . , 
then the following expansion can be used: 

m oo 



n^i^-s^+i/)/^)' 


n; =1 r(i- aj - + ^-(6 ft + i/)/s ft )" 


(-1)" 


nj =m+ ir(i-^ + J B,(^ + ^)/^ 


n?=„ + ir(%-A j (&, + i ,)/ J B,)" 



X- 



v\B h 



(C2) 



(C3) 



where a prime in the product means j 7^ /i. This series converges for all z 7^ if fi > 0, and 
for < \z\ < 1//3, if = 0, where \i and (3 are defined as 

<? v 



v q 

Ai TT n -6i 



/' » »• ; n v 



(C4) 



3=1 3=1 
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II) If the poles of 



JJra-a,--^-) 



are simple: 

A h (l - aj + Aj(l -a h + A), j ^ h, j,h = \,...,n, A, v = 0, 1, 2, . . . , 
then the following expansion can be used: 

n oo 

W z )=EE 

ft=l i/=o 

n;=i r(i - s - - ^ + [n; =1 + - ah + ^m,/ 



n? =m+ i r(l - - 5,(1 - a, + [n? =n+1 r(a, + Ml -a h + v)/A h ) 

( — \y(±}()--ah+v)/A h 



v\A h 



(C5) 



(C6) 



(C7) 



where a prime in the product means j 7^ /i. This series converges for all z 7^ if /x < 0, and 
for |z| > 1//3, if /i = 0, where \i and /3 are defined as above. 



Appendix D: The Mittag-Leffler Function 

To evaluate the Mittag-Leffler function we concentrate on the integral [Eq. (25)] 

e~ xt x a ~ 1 dx 



F a (a;t) 



a sin an f c 
v Jo 



(Dl) 



x 2a — 2a cos an x a + a 2 
Since the Mittag-Leffler function represents behavior between a power law and an exponen- 
tial, it plays a crucial role in both scientific and engineering applications. In this regard, it 
is imperative that we have a closed expression for this integral, which in previous works [18, 
22-24] was assumed to be the purely decaying part of the Mittag-Leffler function with an 
imaginary argument, hence neglected in subsequent calculation with serious consequences. 
Since the standard methods do not allow us to evaluate this integral, we are going to use of 
the language of if— functions. 

We start by noting that the above integral requires the evaluation of the Laplace transform 



£ 



x 



a-l 



x 2a — 2a cos an x a + a 2 



We first express the function 



1 



x 2a — 2a cos an x a + a 2 ' 



(D2) 



(D3) 
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in terms of H— functions. For this, we evaluate its Mellin transform: 

poo 

M{f(x)} = f(s)= f{x)x s - l dx. (D4) 
Jo 

Using Equation (13.126) in Bayin [6], we evaluate the above integral to obtain the Mellin 
transform 

M 



x 2 — 2a cos an x + a 2 
which after using the property 



7r(— a) s 2 sina(s — \)n 



sm7rs sin cm 



M{f(xP)} = ±f(p, (3>0, 



allows us to write 



M {/(*)} = 



n(—a) s / a sin(s — a)n 
aa 2 sin an sin — 

a 



We now use the relation [6] 



7T 



sm7T2; 



to write 



h(s)=M{f(x)} 



T(z)T(l- Z y 
n (-a) s / a T(s/a)T(l - s/a) 



a<7 2 sina:7r T(a — s)T(l — a — s) 
The inverse Mellin transform, which is defined as 

f 7+ioo 



9( x ) = it-. I g(s)x- s ds, 



where g(s) is the Mellin transform of g(x), allows us to write 

1 /IT 



x 2a — 2a cos an x a + a 2 



aa" sin an 



x 



1 p +to ° (-a) s / a T(s/a)T(l - s/a) 
2n~i 



'-/-too T(a- s)T(l-a- s 
Comparing with the definition of the H— function in Appendix (B): 

rrm,n ( ,(a 1 ,A 1 ),...,(a p ,A P )\ _ J_ f h( s ) z 'Sj 

P'i ^I(6i,Bi),...,(6„b,) > ; - 2m J c { ) ' 
we read the values of the parameters as 

m = 1, n — 1, p — 2, q = 2, 

a\ = 0, A\ = 1/a, a,2 = a, A 2 = —1, 

6i = 0, B x = l/a, b 2 = a, B 2 = -1, 



x s <is 



(D5) 

(D6) 

(D7) 
(D8) 
(D9) 

(DID) 



(Dll) 



(D12) 



(D13) 
(D14) 
(D15) 
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thus obtaining 



x 2a — 2a cos air x a + 



a z \ daysman/ 



(0,l/a),(a ,-l) N 
(0,l/a),(a,-l). 



(D16) 



To evaluate F Q (<r; t), and finally the time dependence T(t) of the wave function, we need 
the Laplace transform 



£ 



x 



a-l 



— ) 

sin air/ 



x 2a — 2a cos air x a + a 2 ) \ < \ a- sin i 

x 



x 



'2.2 1 (.^l/a 

Using the following expression of the Laplace transform 

£ {^XT ( 

with the replacements 

p ->■ a, a ->■ l/(-o-) 1/a , a -> 1, s -> t, 

m — 1, n = 1, p = 2, q = 2, 

ai = 0, = 1/a; a2 = a, A 2 = —1, 

6i = 0, Si = 1/a; 6 2 = a, 5 2 = -1, 

we obtain the needed transform as 

„a-l 



(0,l/a),(a-l) N 
(0,l/a),(a-l) / 



.X" 



7T 



:r 2a — 2(7 cos a7r x a + a 2 j aa 2 sin «7T t a 
where we have used the relation [25] 



Trm,n ( \{ a p>A p )\ TTn,rn I ^ 

U p,q \ Z Hb q ,B q ) ) ~ U q,P I 



Combining these results, we finally obtain the desired integral: 

e~ xt x a ~ 1 dx 



„ , s ai a sin a:7r /" c 

F a (<r;t) = / 

n Jo 



x 2a — 2a cos an x a + a 2 ' 



in closed form as 



aat a ' 3 V 7 l(a,l),(0,l/a),(l-a 



(D17) 



)}=s- p H;r i : l {as-%J^), (D18) 



(D19) 

(D20) 
(D21) 
(D22) 



(D23) 



(D24) 



(D25) 



(D26) 
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The time dependent part of the wave function [Eq. (24)] can now be written as 



T(t) = T(0) 



■ .,1/c*. 



a 



1 -t%\(< _\l/a .l(0,l/a),(l-a-l) 



^2,3 U °) r |( Q ,l),( ,l/a),(l- 



•a-1) 



(D27) 



We now use Equation (C3) to write F a (a; t) in computable form as 

FaM = _if ( _ ir ..^,£(zpi±l)(Mg _ £ %!!)_. (D28) 



u=0 v y ^=0 

Since /i — 1 > 0, the above series converges for all |z| 7^ 0. In the first series, we concentrate 
on the expression 

/ = ( iy(i+j:) r (~«) r ( 1+ «) . 



We extend the formula 



rf-n) 



= (-1) 



r(-z/)z/! 

AT! 



N—n ' 



T(-N) nl 
where N and n are positive integers, to non integer arguments as 

r(-g) _ r nQ -, r(Q + i) 
r(-Q) 1 J r( g + i)' 



and write 



r (-^_ T( _i r - r ^ + 1 ) 



(D29) 



(D30) 



(D31) 



(D32) 



r(-iz) lv ' r(| + i)- 

We have inserted the =F sign to display the ambiguity in the gamma function for the negative 
integer values of its argument, where it diverges as =Foo. Thus, / is evaluated as 



j = = F (-l)^)(-i)5-^ r ( y + 1 ) r ( 1 + 9 



Substituting this into F a (a;t) [Eq. (D28], we obtain 



(D33) 
(D34) 



F a {a,t) = ±-^— } l^ T{1 + ai/) 

u=0 v > 



a * — ' v\ 

u=0 



a * — ' v\ 

u=0 



u=0 



which is nothing but 



F a (a;t) = ±- 



a 



E a (\ n i a t a ). 



(D35) 
(D36) 

(D37) 
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Substituting this into Equation (24), we write the time dependent part of the wave function 
as 



T(t) = T(0) 



., l/a. 
iXJ t 



n i\J t 



+ E a (X n i a t a ) 



(D38) 



a a 

To be consistent with the robust result in Equation (21), we pick the minus sign. Thus the 
time dependent part of the wave function is again obtained as 



T(t)=E a (X ni a t a ), 
where without loss of any generality we have set T(0) = 1. 



(D39) 
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